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Abstract. We consider rectangular random matrices of size p x n belonging to 
the real Wishart-Laguerre ensemble also known as the chiral Gaussian orthogo¬ 
nal ensemble. This ensemble appears in many applications like QCD, mesoscopic 
physics, and time series analysis. We are particularly interested in the distribution 
of the smallest non-zero eigenvalue and the gap probability to find no eigenvalue 
in an interval [0, £]. While for odd topology v = n — p explicit closed results are 
known for finite and infinite matrix size, for even v > 2 only recursive expres¬ 
sions in p are available. The smallest eigenvalue distribution as well as the gap 
probability for general even v is equivalent to expectation values of characteristic 
polynomials raised to a half-integer power. The computation of such averages is 
done via a combination of skew-orthogonal polynomials and bosonisation meth¬ 
ods. The results are given in terms of Pfaffian determinants both at finite p and 
in the hard edge scaling limit (p —»• oo and v fixed) for an arbitrary even topol¬ 
ogy v. Numerical simulations for the correlated Wishart ensemble illustrate the 
universality of our results in this particular limit. These simulations point to a 
validity of the hard edge scaling limit beyond the invariant case. 
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1. Introduction 

The ensemble of real rectangular pxn matrices W with independent Gaussian entries 
is the oldest example of random matrix theory, introduced by Wishart in the context 
of multivariate statistics p~{. Since then more general ensembles built of complex 
and quaternionic matrix elements have found a wide area of applications, ranging 
from physics and mathematics to biology and engineering, e.g. see Refs. |2 [3] for 
reviews and references. These ensembles got several names and are widely known as 
Wishart ensembles because of its inventor, Laguerre ensembles because of its relation 
to the Laguerre polynomials, and chiral Gaussian orthogonal, unitary, or symplectic 
ensembles hinting to their transformation properties. They are concerned with the 
singular value statistics of W while in the case of the complex eigenvalue statistics the 
name Ginibre ensemble is more common. 

Despite the fact that the real ensembles are more versatile than their complex 
and especially their quaternion counterpart, those ensembles are at the same time 
technically challenging. This is particularly true when correlations among the matrix 
elements are introduced. The introduction of correlations can be done in three ways, 
either by a correlation matrix resulting in the correlated Wishart ensemble as it is the 
case in the analysis of real time series [E IE IE 0 IE IE 0U1 [HI HE HE HE HE HE HE HE • by 
extending to non-Gaussian probability distributions of the matrix elements Hi HE HE 
or by adding constraints [IE HE HE • 

When analysing the spectral statistics of the positive definite matrix WW T 
one has to distinguish between two kinds of correlations. The first kind involves 
density correlation functions which can be considered on a global or local scale. The 
second kind comprises correlations only involving individual eigenvalues such as their 
distribution or spacing. These correlations are by definition local objects since on a 
large scale individual eigenvalues with a global separation are usually screened and 
thus uncorrelated. An efficient tool to compute the second kind of correlations are 
gap probabilities, meaning that a certain interval is void of eigenvalues. 

To emphasize the importance of the distribution of individual eigenvalues we 
summarise a few applications: The condition number of a matrix A is the ratio of the 
root of the largest over the smallest non-zero eigenvalues of AA T , which was analysed 
in a random matrix setting in |25j . The smallest eigenvalues are responsible for 
chiral symmetry breaking in quantum field theory where real matrices correspond to 
Quantum Chromodynamics with two colours [26]. Those eigenvalues are very sensitive 
for fitting lattice QCD data to random matrix results, see [27l [28l[80l [32l l83l[84l l35l 136] 
and in particular EE EE for our symmetry class, for the discussion of the importance 
of individual eigenvalues in fitting lattice data. In this context, the number of 
zero-eigenvalues n — p = v > 0 corresponds to the gauge field topology, see m 
for a review on this topic. A further application of these quantities can be found 
in studying topological insulators, see [38] for a recent review. In multivariate 
statistics the smallest eigenvalue plays an important role in high dimensional inference 

EE HE EE HE- 

In the uncorrelated Gaussian case all density correlation functions are known 
most explicitly. For finite p and arbitrary v the fc-point correlation functions of 
the real Wishart-Laguerre ensemble are given by a Pfaffian determinant of a kernel 
involving skew-orthogonal polynomials [261 143] . These are expressed in terms of 
Laguerre polynomials. In the limit p —> oo the local kernels are universal and are 
given by the corresponding Bessel-, sine- or Airy-kernel, for the hard-edge, bulk or 
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soft-edge scaling limit, see [33] for the corresponding expressions and references. We 
are particularly interested in the Bessel kernel in the microscopic origin limit (hard 
edge). The universality of this kernel was shown for non-Gaussian ensembles in [45] 
and for the some kinds of correlated Wishart ensembles in mm- 

Because the individual, e.g. smallest eigenvalue distribution, can be in principle 
expressed through a Fredholm Pfaffian [46] (see also m) of the very same kernel this 
universality is inherited by the individual eigenvalue distributions. Due to this fact we 
can restrict ourselves to the Gaussian case. Apart from this relation the universality of 
the smallest eigenvalue at the soft edge has been proved explicitly m■ Furthermore, 
it was shown in [23] that the known expressions for the smallest eigenvalue in the 
Gaussian ensemble and the one with a fixed trace constraint [22] agree hinting to 
stronger universal with respect to non-differentiable deformation. 

What is know more explicitly about the distribution of the smallest eigenvalue? 
Closed expressions in the quadratic case v = 0 were derived in [33 (c.f. ;4Sj). In [25 
an exact recursive scheme in the matrix size p was set up, leading to closed results 
for arbitrary p for v = 0,1, 2, 3 only. A distinct structure including both polynomials 
and Tricomi’s confluent hypergeometric functions was observed to hold for v even and 
odd, respectively. Exact results for finite p and in the hard-edge scaling limit were 
derived in [50] for all odd v fixed and extended to the fc-th smallest eigenvalue in [30]. 
Both results are represented as a Pfaffian determinant, a structure that was unnoticed 
in the recursion of [25]. Simple approximations following the idea of a Wigner surmise 
were tested in [5TI . An efficient numerical algorithm was used to directly compute the 
Fredholm Pfaffian expression in [36], and most recently an extension to the correlated 
Wishart case was presented in [14] . 

The goal of the present article is the derivation of a Pfaffian representation for the 
case of even v at finite and infinite p announced in [5 2\ ■ Thus we aim at completing the 
picture of the smallest eigenvalue distribution in the real Wishart-Laguerre ensemble. 
The idea in [50] [30] that we will use is to represent the smallest eigenvalue and the gap 
probability as an expectation value of powers of characteristic polynomials. For odd 
v half integer powers appear, which is the technical problem we have to solve. The 
tool we apply are skew-orthogonal polynomials with a non-standard weight function 
containing a square root. The determination of these polynomials is then based on 
the method of Grassmann variables and bosonisation [53} [54], [55], a particular case of 
the supersymmetry method, see [56] HZ] and references therein. 

The relevance of such expectation values of characteristic polynomials including 
half integer powers has been advocated independently and solved in a few special 
cases in [58] . motivated mainly from applications to Quantum Chaos. There, the 
Gaussian orthogonal ensemble is considered. We present results for the chiral Gaussian 
orthogonal ensemble (chGOE) for characteristic polynomials raised to an arbitrary 
half-integer power. 

The outline of the work is as follows: In section [2] we formulate the problem 
in terms of expectation values of characteristic polynomials. These can be computed 
by introducing non-standard polynomials that are skew-orthogonal with respect to a 
weight function containing a square root, see section [3] When expressed in terms of 
these polynomials both the gap probability and the smallest eigenvalue distribution 
exhibit a Pfaffian structure. The building blocks appearing in these expressions, the 
partition function, the polynomials and their kernel are computed in section |3] Here 
we also summarise our exact results for arbitrary p and odd u. In section[5]we take the 
microscopic origin limit (jp — > oo and v fixed and even). Furthermore, we illustrate 
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our findings by numerical simulations, including the correlated Wishart case which 
follows the same universal predictions, see section [G] Our conclusions and discussion 
of open problems are presented in section [TJ 


2. Formulation of the Problem 


We consider the singular value statistics close to the origin of the real Wishart-Laguerre 
ensemble. To this end, we take W to be a rectangular matrix of size p x n with 
“rectangularity” v = n — p > 0 and real entries Wij € R. In QCD the “rectangularity” 
v is identified with the index of the Dirac-operator and, thus, with the topological 
charge of the gauge field configuration, c.f. [37] . The entries of W are drawn from a 
Gaussian distribution with row-wise correlations [5], 

P(W\C) ~ exp i tr WW T C~ 1 ^j . (2.1) 

Sometimes also doubly correlated Wishart random matrices are considered to model 
spatio-temporal correlations, e.g. see ® m nn na- The singular value statistics of 
W are completely determined by the eigenvalue statistics of the matrix WW' which 
is called the Wishart correlation matrix. The measure on the space of rectangular 
matrices d [W] is the flat measure, the product of all independent differentials. In 
all of our analytic computations we take C = l p , i.e. we consider the uncorrelated 
Wishart model. Only in section[G]we argue that generically also the smallest eigenvalue 
of the correlated model follows the universal distribution derived in section [5] 

In order to consider the statistics of the eigenvalues of the Wishart matrix we 
diagonalize WW T = OXO T , where X = diagpi,... ,x p ) > 0 and O € O (p) is an 
orthogonal matrix. This leads to the following normalised joint probability distribution 
function (jpdf) of the eigenvalues, e.g. see [55] 00] 

P(X) = J- | A p pO| n 4"~ 1)/2 exp i-Xi/2) . (2.2) 

1=1 

The constant Z v ^ v is the partition function and, hence, the inverse of the normalisation 
constant. It is a Selberg integral [59] and explicitly reads 


4/ iA„mi=ff r[(,+ 3 )/ 2 (p l / (j 2 +1/+1)/21 

*-ln J—0 


(2.3) 


for arbitrary v > — 1. The term A P (X) = II — x i) denotes the Vandermonde 
determinant. The variables Xi coincide with the squares of the singular values of W 
and typically describe the low lying eigenvalues Ai^fxl of the QCD-Dirac operator 

mm 

For the smallest non-zero eigenvalues two kinds of large-p limits have to be 
distinguished. If we take n and p to infinity while c = p/n, c £ (0,1], is kept fixed, the 
macroscopic density of Marchenko-Pastur detaches from the origin. The level density 
vanishes with a square root behaviour at both endpoints which are called soft edge. In 
this situation both, the largest mm and the smallest eigenvalue, are Tracy-Widom 
distributed SH- 

In contrast to this soft edge scaling we can fix the index v = n — p when taking n 
and p to infinity. Then the macroscopic density behaves as an inverse square root at the 
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origin, also known as the hard edge. The corresponding scaling is called microscopic 
origin limit. This limit will be considered in section [5] The microscopic level density 
HUES and all fc-point correlation functions are given by the Bessel-kernel [33] which 
is universal [331. In principle all individual eigenvalue distributions including the 
smallest follow from these density correlations for arbitrary fixed v. They are given by 
the Fredholm Pfaffian of the Bessel kernel [46i. In particular the limiting distribution 
of the smallest eigenvalue, P v (x), follows the microscopic spectral density, p v {x), for 
small x, in particular P^{x) ss p v (x). The level density of the matrix WW^ is [61] 


Pi'(u) — ^ (My/u) Ji/ — 1 {p/u) (\Zu)) + —^'(N^) 


VU \ 

1 - J dsJ„(s) (2.4) 


with J„ the Bessel function of the first kind. 

At finite p and v the distribution of the smallest eigenvalue, P p>!/ (f), can be 
derived via the gap probability E PtV (t). This gap probability is the probability to find 
no eigenvalue in the interval [0, t\. Starting from the jpdf (12.211 we immediately have 


E P ,v(t) 


p 

-J— n [ dx i x t~ 1)/2 e- Xi/2 lAppOl 
^ *= 


(2.5) 


o-pt /2 


Zr) 


p p 

[] J d Xi (xi + |A P (V)|. ( 2 . 6 ) 


In the second line we have shifted Xi —► Xi + t for all i = 1 leaving the 

Vandermonde determinant invariant. 

Let us define the expectation value for any integrable quantity f(X) that only 
depends on the eigenvalues, 


1 p r 

</P Q) p , v = y- n / d ^r 1)/2 e- Xi/2 f(X) |AppOI • (2.7) 

J 

i—i 0 

Then the gap probability (EH) can be written as an expectation value of a 
characteristic polynomial to a certain power with respect to a different partition 
function with fixed index v = 1 , 

E Pt „(t) = e -? t / 2 ^(de&- 1 V 2 (X + tl p )) . (2.8) 

For v = 2k + 1 odd the determinant is raised to an integer power 7 = (z^ — l)/2 = 2k. 
Such expectation values have been computed for finite dimension [50] as well as in 
the microscopic origin limit [30]. For v = 2k even the expectation value is taken of 
a half integer power of a determinant. This average was up to now an open problem 
and is computed in the present work. The special cases with v = 0 [49] and v = 2 
[23] were computed in the microscopic limit using a different route, where the latter 
result is based on the recursive construction of Ref. [25] that only yields closed form 
expressions for finite p at v = 0 and v = 2 even. 

Next we turn to the distribution of the smallest eigenvalue. This distribution 
follows from the gap probability by setting the smallest eigenvalue equal to t , say 
X\ = t. Hence it is the first derivative in t of the gap probability (12.51) . 

t 

P p ,„(t) = ~^ E pAt) ** E PtV (t) = 1 - J dt'P p ,„(t'), (2.9) 

0 
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where the second relation fixes the normalisation. Rather than first computing 
and then differentiating with respect to t we can directly start with P p v (t), 


P p ,v{t)=——& 1)/2 e 4/2 / d Xixf 1)/2 e Xi/2 {xi - £)|A p _i(xj> 2 )| 

^ i=2 J 

p °? 

= / dxiXi(xi +t) l ' u ~ 1)/2 e~ Xi/2 \A p - 1 (xj> 2 )\, 

^P,v ] — 9 j 


Ep, v {t) 


( 2 . 10 ) 

( 2 . 11 ) 


where A p _i(xj> 2 ) does not contain x\. We recall that the smallest eigenvalue is 
chosen as x\ = t. Therefore the absolute value of these terms can be dropped. In the 
second line we have again shifted Xj —> Xi + t for all i = 2,... ,p. Consequently also 
Pp,u(t) can be written as an expectation value, 


Pr 


At) = p^- 1 )/ 2 e-^ 2 %^ (det ^ 1 )/ 2 (X + tl p ^)) 

Z p , v \ / p-1,3 


( 2 . 12 ) 


Here the expectation value is with respect to a partition function of p — 1 eigenvalues 
with v = 3 fixed, accounting for the extra factor Xi in Eq. (12.111) . 

Once again for v = 2k + 1 the expectation value of determinants to an integer 
power k is known. Our task reduces to the computation the half integer case v = 2k. 
Summarising the problem we consider the following type of expectation values 


/def 1/2 (X + fl p ) J|det(X + tjl p )\ . (2.13) 

\ 1=1 / 

Such problems were advocated independently in Ref. [55] for the Gaussian orthogonal 
ensemble also including more than one square root in the denominator. While special 
cases have been computed in Ref. [58] . having different applications in mind, we 
determine Eq. (12.131) . leading us to the gap probability and the smallest eigenvalue 
distribution at finite matrix dimension and in the microscopic limit. 

As a final remark the gap probability (G3D has been studied in more detail 
for even rectangularity v = 2k in M, where a dual supermatrix model was 
found. Although the supermatrix model is invariant under the action of a particular 
supergroup, it was not solved in [14], because of non-trivial subtleties related to the 
necessary diagonalisation. We circumvent these subtleties by combining the method 
of Grassmann variables with the theory of orthogonal polynomials. 


3. Pfafflan Structure and Non-Standard Skew-Orthogonal Polynomials 

We tackle the expectation value (12.131) by including the unwanted inverse half-integer 
power of the characteristic polynomial into the weight function. The remaining 
integer powers can be expressed in terms of skew orthogonal polynomials (SOP) using 
standard techniques. The difficulty is thus shifted into finding the SOP with respect 
to the t-dependent weight 

wAx,t) = ^==e- vx/2 - (3-1) 

y/X + t 

The exponent 7 is equal to 7 = 0 for the gap probability and to 7 = 1 for the smallest 
eigenvalue. Thus we calculate the SOP in a unifying way for both quantities. The 
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auxiliary parameter y is set to unity unless otherwise stated. It is needed to generate 
the polynomials of an odd order. 

We seek monic, parameter dependent polynomials R^\y,t) = yl + ... that are 
skew-orthogonal with respect to the anti-symmetric product 

oo oo 

y — 3 ? 

i _ i w^(x, t)w 1 (y, t)f(x)g(y) 

0 0 



oo y 

= I dy I dxWr y( x ’ t ') w 'y(y’ t ') (f( x )s(y) ~ f(y)3( x )) ( 3 - 2 ) 

0 0 

for two arbitrary integrable functions / and g. In particular the polynomials have to 
fulfill the relations 




( 7 ) T}{~t)\ _ 


R\i’R. 


= (r: 


?(t) 


,RilL) =0 


.... x x (3-3) 

t J \ J / t \ J It 

for p = 2L even. The parameter dependent constants r,- 7 ^ (t) are their normalization 
constants. In the case p = 2L + 1 odd the polynomials have to satisfy an additional 
condition [SSI 


= r ™ (t)s 
dxR(x, t)wj (t; x) = 5i t 2K, 


5 ( 7 ) \ 


„(7 )/ 


Jij j 


RiV,R h) \ - 




_ / r?0) BO) 

2i / t ~ \ n 2j+V n 2i+l / t 


= 0 , 


(3.4) 


which fixes the normalization of the polynomial of highest order R^-(x, t) with K > L. 

Note that the normalization constants rj 7 '* (t) are the same as for the case p = 2L. 
The reason is the following relation between the two different kinds of the polynomials 


R { ?\y,t) = R^\yR) - 


fo°° d xw 7 ( x , f)i?.j 7) (x, t) n(7) 


/ 0 °° d xw 1 (x, t)R^(x, t) 




(3.5) 


for j < 2K and R^{y,t) = 


- pM 

2 K 


R K or(y,t)/ f™ dxw-,(x,t)R^(x,t). 
?(7 )/ 


('>') (m- A Therefore 
We emphasize that the 


we concentrate on the polynomials R ' f {y, t) only. 

additional condition / 0 °° dxR l f l \x, f)w; 7 (f; x) = 8i^K can be also replaced by other 
conditions. For example in the framework of [62] the condition would read 
J 0 °° dxR^\x, t)uij(t-, x) = Si^o for all i > 0 which has other advantages in the 
calculation. 

We define a new partition function 


p 4r 

Z p ^(t) = Y[ J dxiWry(xi]t) | A p (X)| 

/»00 

< dxw 7 (x, t)R^{x, t), 


LP/2J-1 

= P- n r 5 7) 

3=0 


ry it) 


7 \ 1 J 2L ' 

1 , 


p = 2L + l, 
p = 2L 


= Z Pt 27+1 (def 1/2 (X + tl p )) , (3.6) 

where the second line is a general relation between the normalization constants of the 
SOP and the partition function [59]■ The floor function |_p/2j yields the largest integer 
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smaller than or equal to p/2. Employing this partition function we introduce a new 
parameter dependent expectation value of an integrable observable /, 

P 00 

</PO>p, 7 = ~"7TT n [ d XiW 7 ( Xi -,t) f(X ) | A P (X)|. (3.7) 

i-i Q 

The parametric dependence on t is indicated through the superscript. 

In this framework the gap probability (12.81) and the distribution of the smallest 
eigenvalue ( 12 . 121 ) (both for v = 2k) read 

E P Mt) = (det fc (X + tip))' , (3.8) 

PpMt) = pt ^-^e-^ 2 Zp - lAt) Uet k (X + tip.,))' . (3.9) 

X/p,2fc ' ' P~ 1,1 

It is worth emphasizing that now only integer powers appear in both expressions. 
Thus, we can apply the results from the literature for general weight functions, c.f. 
mmm- For k = 2m even, we obtain the following Pfaffian expression with a 
2m x 2m dimensional kernel, /C p +fe(«: a , n b ), 

I U \ t 

( P + n\ k 7^]t\ A \ [£p+k(«a, «&,*)] , (3-10) 

(P + k)\Z Pi7 (t) A k (n) 


7 det(X - K a l p ) 


\ a—1 


P.7 


while for k = 2m + 1 we have 
( det(X - k 0 1 p )\ 


P^Zp+k+l,~/{t) 


(p + k + l)\Z pn {t) Afc(«) 


(3.11) 


x P f l< 0 


fcp+k+l (^a? ^b) P p-t-fc+1 (^a; Z) 


" b<k -Tp+k+l(Ka,t) 0 

The kernels inside the Pfaffians are given by 


B%0 + 1 (/C 0 ,i) («&,*) - -R^+l («6, *) ( K a,t) 


ICl(K a ,K b ,t) = < 


E 

i=o 


2 j 


'2j+l 


2 j 


rj 7 ) (t) 


Z G 2N, 


(( 3)/2 _g(7^ ( Kat f) («{,, f) - ^2]+! *) ( K a, Z) 


E 

J=0 


r5 7 ) (t) 


ZG 2 N+ 1 . 

(3.12) 

The case of k = 2m + 1 odd is obtained here from the case k = 2m + 2 even 
by introducing an additional determinant in the average depending on the dummy 
variable « 2 m+ 2 - This variable is sent to infinity such that the additional row and 
column in Eq. (13.111) reads 

Xli^Ka^ ^2m-f-2, t) 


Ri{n a ,t) = - lim 

«2m+2—>00 


K 1 - 1 

ft 2m+2 


(3.13) 


This limit is independent of l being even or odd. 

Once we have determined the SOP and their kernel we have to take the limit, 
n a —> —t for all n a - This yields derivatives of the polynomials R ■ 7 - ) (n b ,t) because 
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of l’Hospital’s rule. Moreover we could also just take part of the K a —> —t which 
is needed to calculate the distributions of the second smallest eigenvalue, the third 
smallest eigenvalue etc., see Ref. [503■ The distributions of the smallest eigenvalue for 
QCD with dynamical quarks can be found in this way, too, cf. Refs. [30:. 

In the next section we explicitly compute the SOP. For this computation it is 
helpful to understand also these polynomials as expectation values [65; ;66j, 


R 2 j\v^) = (det {yl 2j - , (3-14) 

= (det (yl 2 j+i — X) (y + Cj(t) + TrX))*^ +1 ^ 

= (y + c 'j (*) - 2 ^) r( S (y> *) ^ 

= (y + Cj{t)-2y-^--2t^j R { ^\y,t). (3.15) 


These two last relations also hold in a much more general framework where 
w 1 {x,t)w 1 {y,t){y — x)/\y — x\ is replaced by an arbitrary anti-symmetric two-point 
weight g{x,y) = -g(y,x) [67]. 

The odd polynomials are not unique [65] [59] which is reflected in an ambiguous 
constant cj(t). This constant can depend on t and the index j but is independent of y. 
This dependence is the reason why one could absorb the derivative of the normalization 
into dj(t) and rephrase the derivative in r) as a derivative in t and y , yielding a new 
constant in Cj(t). We will stick to the derivative in p here. 

The kernel can be directly expressed as an expectation value, too. Rather than 
computing the individual polynomials (13.141) and (13.151) and performing the sum (13.121) 
one can consider the average 




l(l-l)Zj- 2 „(t) 


(k q - Kb) (det (X - Kali-2) det (X - K b li - 2 ))\_ 2 l , 


(3.16) 


e.g see [65] in the hermitian limit or [63] in the general framework of anti-symmetric 
two-point weights. This representation is also useful when proving that the large-p 
limit for even and odd p yields the same answer. The additional row and column in 
eq. (13. lip is then 

t) = l{l ~ l z l Z ^ l{t) (det (X - K a h. 2 ))l^ . (3.17) 

We emphasize that for l even this function is equal to the polynomial i? z _ 2 ( K a), up to 
a constant. 


4. Calculation of the Finite p Results 

We start our calculation by considering the expectation value 

= (H det(X - K a li)\ . (4.1) 

\ &= 1 / 

This quantity is a polynomial in the variables K a ■ The highest power in these variables 
determines its normalization, Ii(k) = (— l) lk K[ ■ ■ ■ K l k + .. ., such that we can omit the 
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normalization constants in the intermediate steps of our calculation. The overall 
constant can be fixed at the end of the calculation. 

In the first step we rewrite Eq. 63 } as an integral over a rectangular real matrix 
W of dimension l x (l + 2y + 1), where in this and the next subsection 2y + 1 can be 
any integer, 


r. k 

Ii(k) oc J d[W] det(« a lz - WW T ) 


exp [~trWW T /2 S j 


OC 


d[W] n det(K a l; — WW 1 


•y/idet ^WW T +1 

exp f— tr WW T /2) 

1 V = 

J det (w T W + il J+27 +i) 


(4.2) 


We emphasize that the normalization constants in each of the steps may dependent on 
t. In the second line of Eq. (14.21) we replaced WW T —> W T W because of the relation 
between the determinant and the trace, i.e. In det A = trlnA, and the invariance 
of the trace under circular permutations. The switching of the order of W and W T 
allows us to avoid the Efetov-Wegner terms appearing in the superspace dual to this 
average, see for details Refs. 5Bj [TOl 17} 1551 . 

We introduce a real (l + 2y + 1 (-dimensional vector v to rewrite the single 
determinant in the denominator as a Gaussian integral. Additionally we express the 
product of determinants in the numerator as a Gaussian integral over a rectangular 
matrix V of dimension l x 2k whose entries are independent Grassmann variables 
(anti-commuting variables). For an introduction in supersymmetry we refer to 
Refs. |70i [56, El] and for the supersymmetry method with general weight to 
Refs. [71115511541 R)51 155] , Particularly the bosonisation is described in Refs. [5511511 1551 . 
The matrix V satisfies the following symmetry under complex conjugation and under 
Hermitian conjugation, 


V* = V 


0 

-1 k 


1 k 

0 


and (P 1 ') 1 ' = — V, respectively. 


(4.3) 


Then, the average reads 


1 r 


Ii(k) oc j d [W] j d[V] / d[u]exp(— ^ tr WW T + tr WW T VV^ + tr Wvv T W T 


x exp 



[tr V^Vk + tv T v] 


(4.4) 


with k = diag(fiq,..., «*;) <S> I 2 - Because of the symmetry (|4.3I) the dyadic matrix VV't 
behaves like a real symmetric matrix such that we can integrate over the matrix W 
without symmetrizing the other terms. This integration yields 

Ii(k) oc J d[E] J d[v]det ” 1/2 (l z 0 1 (+ 2 7 +i + W f 0 1 (+27+ i + 1; 0 vv T ) 

x exp i [tr V^Vk + tv T v\^j 

oc J d[V] J d[v]det -1 / 2 ([1 + v T v\±i + VV^ 

x det - ^" 1 " 27 ^ 2 (lj + VV *) exp - [tr V^Vk + 
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a J d[H] J d[u][l + v T v]-^ l+ 2 k ^ 2 det 1/2 ([1 + v T v]l 2k + V^V) 

x det^ +27 ^ 2 (I 2 k + V^V) exp ^ [tr V^Vk + tv T v]\ . 


(4.5) 


Again we have used the relation between the determinant and the trace and 
the invariance of the trace under circular permutations. However we have to 
remind ourselves that anti-commuting variables are involved such that tr(V’V"r) m = 
— tr(V'V) m for any m € N. This explains the change from negative to positive powers 
of the determinant in VV^. 

In the last step we can choose between two approaches, the generalized Hubbard- 
Stratonovich transformation no ng and the superbosonization formula [ 551 1541 . 
Both approaches are equivalent [55]. We choose the superbosonization formula 
since it directly leads to a compact expression. Since no supermatrices comprising 
both, bosonic and fermionic, blocks are involved the superbosonization reduces to 
bosonisation, only. This means that the norm v T v is replaced by a positive variable r 
(the square of the radial part of an ordinary real vector) and the dyadic matrix V'V 
is replaced by a self-dual, unitary matrix, 


U = 


0 

-1 k ' 

U T 

0 

l fc ' 


0 


0 


G U(2 k), 


(4.6) 


because V^V is itself self-dual. The set of matrices defined via Eq. m is the 
circular symplectic ensemble first studied by Dyson [72]. This set is the coset 
CSE(2fc) = U(2/c)/USp(2fc) and has a uniquely induced Haar measure dy(U) from 
the unique, normalized Haar measure of the unitary group U(2fc). In particular up 
to a normalization constant it is given by dp(U) oc dfUJ/det^ 1 / 2 !/, with d[f7] the 
product of differentials of all independent matrix entries of U. The superbosonization 
formula yields 

r r°o ( l - |-27-l)/2 

/,(«) = Cf 1 J d M (I/)det-V 2 (/det<‘^>/ 2 (1„ + t/)jf d q 1+r) |, +2t)/2 

CSE(2fc) 

x det 1 / 2 ([1 + r]l 2 fc + C/)exp ^ [tr Un + tr]^ , (4.7) 

with the normalization constant 

r , r°° r-d+27—1)/2 

Ci= / d^{U)det~ l/2 Ue tlU / 2 / dr—-^-e-*^ 2 . (4.8) 

J Jo (1 + ry! 1 

CSE(2fe) 


The constant follows from the asymptotics for k —> 00 . The powers of the additional 
terms det~ f ' 2 C/ and r^ l+2 ' 1 ~ 1 ^ 2 only reflect the nature of the variables from where U 
and r originate. We underline that the half-integer of the determinants do not cause 
any problems since the matrices are Kramers degenerate. Therefore the determinants 
of them are exact squares and the square root is taken such that the result is still a 
polynomial in the matrix entries. 

Starting from expression (1471) we calculate the partition function Z p ^(t) in 
subsection 14.11 the polynomials R^\y,t) and the function iFi(Ka) in subsection 14.21 
and the kernel K.i{n a , Kb) in subsection l4.3l In subsection l4.4l we collect everything and 
give explicit expressions for the gap probability and the distribution of the smallest 
eigenvalue. 
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4-1. Normalization Constants 

The partition function Z prj (t) is equal to the case l —> p and k —> 0 in the integral (14.71) . 
Therefore we have no integral over a circular unitary ensemble and only the integral 
over r remains, i.e. 

( p +2 7 —1)/2 


,~tr /2 


(4.9) 


z ( t )= _^ r dr ^z 

mU 2(P+^+ 1 )/2r[(p + 2 7 + l)/2] Jo (1 + r)p/ 2 

The constant is fixed by the asymptotic behaviour Z pn (t) = Z p ^ 1 +it~ p ^ 2 + o(t ~ p / 2 ) 
for t —> oo. The remaining integral is a Tricomi confluent hypergeometric function m, 


i r°° 

\J(a,b,t) = ——- / dzz a ~ 1 (l + z) b ~ a ~ 1 e~ tz , with Rea, Ret > 0 and a, &, t G C. 

i>] Jo 


This hypergeometric function was already found in the work by Edelman 
a crucial ingredient in his recursive formula. The partition function reads 

z pl {t) = 2(p- 1 )(p+27+i)/2 ( HO’ + 3 )/ 2 ] r [(j +2 7 + 2)/2] 

V'=0 r I 3 / 2 ] 

x ^+1/2 V fP± 27+1 ^7 + 3 t\ 

\ 2 ’ 2 ’ 2 J 


(4.10) 
and is 


(4.11) 


(p -1 


_ 2p(p+ 2 7)/2 I J 

+'=° 


T[(j + 3)/2 ]T[(j + 2 7 + 2)/2] 

r[3/2] 


'p 1-27 t 

2 ’ 2 ’2 


The second equality follows from the Kummer identity of Tricomi’s confluent 
hypergeometric function [73], U (a,b,t) = t 1_b U(a + 1 — 6,2 — b, t). Combining (12.31) 
and (13.61) we have thus obtained the first building block for the Paffian structure, 

(4.12) 


(detr I/2 (x + 11,)/ ' = 2-P /2 U f ■ 

which is even valid for any 1 / G No- 

In the particular case of the gap probability the constant is 

(11 j ^(^-14)= 2P/a (n/) u (|.-n 

(4.13) 

We have used the duplication formula of the Gamma function, r[z]T [2 +1/2] = 
2 1 ^ 2z v / 7rT[22:], to simplify the expression. Also the partition function needed for the 
distribution of the smallest eigenvalue the partition function takes a simple form, 

v.,.(<)=2 <p - 4 ' /2 (nil) c /2 u (^.5,1) 


\3= 

P 


= 2 (p-D/ 2 [ J]j! ] U 
U =1 


p — 1 It 
2 ’ ~ 2 ’ 2 


(4.14) 
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The second expression of Tricomi’s confluent hypergeometric function is employed in 
Edelman’s work [25] . 

The normalizations r^\t) of the SOP can be deduced by combining Eqs. (13.61) 
and (Id. 111 ) . 

Z2j+2,y(t) 


^\t) = 


( 2 j+ 2 )( 2 j + l)Z ail 7 (t) 

U (j + 7 + 3/2 ,7 + 3/2, t/2) 


= 2(2j)!T[2i + 2 7 + 2] 
= 2(2j)!T[2j + 2 7 + 2] 


U(j +7 + 1 / 2,7 + 3/2, t/2) 
U(j + 1,-7 +1/2, t/2) 


U (j, -7 + 1/2, t/ 2 ) 

This term becomes important in the sum in the kernel. 


(4.15) 


4-2. Skew-Orthogonal Polynomials 

First we concentrate on the function lFi(K a ,t), see Eq. (13.171) . For this case we set 
l —> l — 2, k —> 1, and n —> n a l 2 . The 2x2 dimensional unitary matrix U only consists 
of one phase e l(p on the diagonal due to its self-duality. Thus we have to calculate the 
double integral 


Pi(n a ,t) oc d p 


Z7T 

/' 


(l + e+) 


Z+27-2 


2 )<p 


> r (Z+2 7 -3)/2 

iT 


tr / 2 (l + r+ e+) . 


(4.16) 

The term (77 + r+ e* v ) is the only coupling between the two integrals and yields a sum 
of two terms, of which each is a product of two functions. The integrals over r are 
equal to Tricomi confluent hypergeometric functions (Id.101) while the integrals over ip 
are modified Laguerre polynomials in monic normalization, 

2 t r 


L^\y) = (-l) Q a! J ^e“ a ^(l + = y a + ... (4.17) 

0 

Then the function appearing in the additional row and column for odd /c, cf. Eq. (13. lip , 
is 


Flfaait) = (- 1 ) 


t l(l - l)Z;_ 2 , 7 (t) f T ( 2 7 ) 


Zl,y(t) 


Likina) 


(4.18) 


n 0^7 + ~ 1 )/ 2 '7 + l/ 2 ^ t/ 2 ] T (27+1) [ , 

( U[ 7 + (Z — l)/2, 7 -(-3/2, t/2] l ~ 3 { a) 

The function was normalized via the known expansion to leading order in n a 

Ti( Ka ,t) = (-l) l l(l-l)Zi-2, 1 (t)K l ~ 2 /Z l , 1 (t) + ... (4.19) 

From the expression (14.181) we can readily read off the polynomials of even order 

(l —» 2 j + 2 ), 

R h) (v t) = L (27) (V) - 2l - U b' + 7+ l/2,Z/2 ] (27+ d 

% (y, ) 2, (y) + 1/2,7 + 3/2,t/2] L «- 1 


) V — 2j 

_ U[j + 7 + 1/2,7 + 1/2, t/2} ^( 2 7 +i) 


(v) 


U[j +7 + 1/2,7 +3/2, t/2] ^ 

2 j + 2y + l U[j +7 + 3/2 ,7 + 3/2, t/2] ( 2 7 ) 


'U[j +7 + 1 / 2,7 + 3/2, t/2] + 


L\ 7 \y), 


(4.20) 
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for any integer j > 1 and R^\y,t) = 1. The second expression is the one presented 
in Ref. [52] and can be found by splitting the term (1 + r + e l(p ) in (1 + e llp ) and r 
instead off e llp and (1 + r). 

The odd polynomials are determined by their relation (13.1511 to the polynomials of 
even order. For this purpose we recall some recurrence relations of the monic Laguerre 
polynomials and the Tricomi confluent hypergeometric functions, 

(y - 2 y§^j L( a\y) = L { ah(y) + + l ) L[ a\y) - + ^L^^y), ( 4 . 21 ) 

U (a, 6, 0 = (1 - b) U (a, b, 0 + (b - a - 1)U (a, 6 - 1,0. (4.22) 

These relations yield 

K%j+i(y, = L 2 j+i(v) + ( 2 7 + 1 +c j (t))L^ 1 \y) - 4j(7 + j)L 2 j-i(j/) 

+ (i)Tg 7+1) (y) + df (t)Lg1+ 1} (2/) + 4 3) (2/) ( 4 - 2 3) 

with the coefficients 


df\t) 


. U[j + 7 + l/2, 7 +l/2,t/2] 
? U[j + 7 + l/2,7+3/2,t/2]’ 


(27 + Cj(t))df ) (t) + (4 1} W) 
-2(2j-l)( 7 + j)4 1) (t). 


4j(i + i) 


U[l + 7 + 1 / 2 ,7 ~~ 1 / 2 , t/ 2 ] 
U[i + 7 + 1 / 2,7 + 3/2, t/2] 


(4.24) 


We have already identified part of the terms in 


df\t) and df\t) 


with the coefficient 


We underline that the ambiguous constant Cj(t) for the odd polynomials, cf. 
Eq. (13. 151) . is not fixed, yet. Thus we are free to choose the coefficient Cj(t) = — 2y— 1 
such that one of the polynomials drops out in (14.2311 . Then the polynomials of odd 
order are a linear combination of only five Laguerre polynomials. This simplifies 
the result presented in Ref. |52]. Nevertheless we emphasize that both results are 
correct due to the various relations satisfied by the modified Laguerre polynomials 
and Triconri’s confluent hypergeometric functions, and the ambiguity in the constant 
Cj(t) which was chosen differently in Ref. [52] compared to the simpler choice here. In 
the microscopic origin limit performed in section [5] we choose another constant Cj (t) 
to simplify the asymptotic result. 

The polynomials R^\y, t), needed for the case p = 2L + 1 odd, can be obtained 
from the polynomials R^\y,t) with the help of the relations (13.51) . Therefore the 
polynomials of even order, R^\y,t), are a linear combination of four Laguerre 


polynomials and the polynomials of odd order, R^ +1 (y,t), can be expressed as a 
sum of six Laguerre polynomials with a suitable choice of the constant Cy(t). 


4-3. Kernel 

The kernel JCi(n a , Kb,t) can be first of all understood as a sum over the SOP, see 
Eq. (13.121) . Plugging the results of subsection 14. 2 1 into this sum we are done. However 
we can also start from the representation (13.16[) and take the general result (14.71) for 
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l —> l — 2, k —> 2, and k —> diag(/t a , K b ) ® 12- Then we find 
K'l (Kq 5 Kfr, l) 

(K a - K b ) 


oc 


j d^([/)det" (i " 2)/2 [/det (i+27 “ 2)/2 (1 4 + U) 


(4.25) 


CSE(4) 


dr- 


(j+2 7 -3)/2 


■det 1 / 2 ([1 + r]l 4 + U) exp ^ [tr Uk + £ r]^ . 


7 0 (1 + r ) d + 2 )/ 2 

To evaluate the integrals we first expand the determinant coupling U and r, i.e. 

det 1/2 ([1 + r]l 4 + U) = (1 + r) 2 + (1 + r) tr 17/2 + det 1/2 17 , (4.26) 

and arrive at 

(Kq 5 Kfo, t) 


OC 


(k q - Kb) 

roo r (l+2 1 -3)/2 e -tr/2 

/ dr- 


(4.27) 


r oo r (l+2j-3)/2 e -tr/2 

dr- 


(1 + r ) 7- 2 )/ 2 i 0 (1 + r )'/ 2 

J dAt(17)def ( '" 2)/2 C/det ( ' +2T “ 2)/2 (1 4 + 17) e“ tr ^ /2 


9 | d 

dK a dKb 


f 


CSE(4) 

°o r (Z+2 7 -3)/2 e -ir/2 

dr- 


J d^(17)def (i " 3)/2 C/det (z+27 “ 2)/2 (l 4 + C/)e- trC/K/2 . 


(X -|- j~)d + 2 )/ 2 

CSE(4) 

The derivatives in and Kf, generate the trace of U. 

Next we diagonalize the matrix U = Vdiag(e* Vl , e l(p2 ) ® I 2 with 

V € USp(4)/USp 2 (2). The normalized measure becomes 

dn(U) = \e iv>1 - e iv>2 fd^d^dp^iV)/^) 2 ) (4.28) 

with d/iHaar(l // ) the normalized Haar measure on USp(4)/USp 2 (2). The integral over 

V is an Itzykson-Zuber integral which is well-known m, 


/USp(4)/USp 2 (2) 


d/in aar (1^) exp 


-tr V&V* 1 k 

2 


(4.29) 


= 6 


( K a - Kb)(e* ¥ ’ 1 - e 1 * 2 ) + 2 _ ivi_ i V2 r / 

e K ° e Kbe + {pi O ^ 2 } 


(k„ - K h ) 3 (e^i - e ^) 3 
We plug this result into the integral over U and have for two arbitrary p, q e N 


J dp(U) det ~ p/2 U det (p+l?)/2 (1 4 + 17) e" tr ^ 


CSE(4) 

2 


r* 27T 


p27T 


di/? 2 r -i(p + 2)(y 1 +<p 2 ) 
27T 


_ d^i 

(k» - Kb) 3 Jo 2?r 
X [(k„ - Kb)(e itpi - e ilp2 ) 2 + 2{e iipi - e iV2 )] e"' 4 '* 
2 1 
(p+ l)!(p + 2)! (Ka - «b) 3 


e 1 ^(1 + e i¥,1 ) p+9 (l + e i¥,2 ) p+ « 




(Ko - Kb) ({p + l)L^ ) (Ka)Tp+ 2 2) (Kfc) - 2(p + 2)(K a )(n b ) 


+ ( p + l ) ip+2 \ K a 


- 2L% 1 \K a )L^\K b ) + 2L^( Ko )L^VM 


-(9-2) 


(9-2)/ 


(<?-!) / 


_d _ d_Y 

([p + 2]!) 2 \_K a - Kb \dK a dKb ). 


Tp+2 2) (K Q )Tp+ 2 2) (Kb). 


(4.30) 























Smallest eigenvalue distribution in the chGOE with even topology 


16 


The second equality results from an expansion of the polynomials in the bracket in 
the phases e llfil and e lv>2 . The third equality is a compact representation by rewriting 
the polynomial in the brackets as derivatives in K a and Kb- We remark that Eq. (14. 3011 
is normalized to the leading order term 2 {K a Kb) p /\p\{p + 2)!]. 

We combine the intermediate result (14.301) . the definition of Tricomi’s confluent 
hypergeometric function (14.101) . and the derivative of monic Laguerre polynomials 
d y lSa\y) = aL^^\y). Then we find for the kernel 


(^a, ZSfr, 


—2,7ft) 

Zl ,7 ft) 


(«a 


Kb) 


1 

K a - K b 




(4.31) 


X 

+ 


\ u(7 +ft —1)A 7 + 1 / 2 , t/2) / 

. U (7 + ft — l)/2, 7 + 3/2, t,/2) \ 
U(7+ft-l)/2, 7 -l/2,t/2) 
u (7 + ft - l)/2 ,7 + 3/2, t/ 2) dK a dK b 


d d 
dK a dKb 


L 


(27-2) 

l 


(k a)L, 


(27-2) 

l 


{K b ). 


Although we have not identified the recurrence relation of the skew-orthogonal 
polynomials we have been able to derive the Christoffel-Darboux formula 
corresponding to the sum (13.121) . It expresses the whole sum as a finite small number 
of terms. Each term is a product of two Laguerre polynomials. The total number of 
these terms is twelve after differentiating and ordering the Laguerre polynomials with 
respect to there index. This is a fixed number of terms which shows that this expression 
of the kernel is ideal to study the large l behaviour including the asymptotic behaviour 
in the bulk and the soft edge scaling limit. Nonetheless we employ the expression as 
a sum over the skew-orthogonal polynomials, see Eq. (13.121) . in the derivation of the 
gap probability and the distribution of the smallest eigenvalue in the microscopic 
limit. The reasons are the additional derivatives we have to perform resulting from 
degeneracy of the variables Kj —► —t, see subsection 14.41 Then the result (14.311) 
becomes quite nasty due to a l/(K a — Kb) term in the differential operator in front of 
the product of the two Laguerre polynomials. 


f.f. Gap Probability and Distribution of the Smallest Eigenvalue 

To obtain the gap probability or the distribution of the smallest eigenvalue itself we 
have to set all variables equal, K\ = ... = Kk = —t. Hence we have to apply l’Hopital’s 
rule in Eqs. (13.101) and (13.111) yielding 

4 (-i)Mfc-i )/yz p+fc , 7 ft) 

{p + k)\Z Ptl {t)Yl k jZoj\ 

X pfl<a,b<fc [^Ki fcp+k{Kl, K 2 , £)|ki=k 2 =—t] 

for k = 2m even and 

t (_i)fc(*-D/2 P \z p+k+ln {t) 


det fc (X + tl p )\ 

/ % 


(4.32) 


(det fc ftA + fl p )^ 


p ’ 7 ip +^+i)^p, 7 ft) rij=o j- 

d^ 1 d b ~ 1 IC p+k+1 {K 1 ,K2,t)\ Kl = K2 =-t d < f~ 1 J r p+fc+l(K,t)| re =-t 
—d b 1 J-pj r k+i{K,t)\ K ——t 


X pf. 


l<a,b<fc 


(4.33) 


0 


for k = 2 m + 1 odd. The additional sign and the product of inverse factorials in front 
of the Pfaffians result from differentiating the Vandermonde determinant Aj.(k), cf. 
Eqs. (13.101) and (13.111) . 
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The prefactor in front of the averages for the kernels K p+k (n\, k 2 , t) and 
polynomials F p+k+ i(n, t), have to be considered together with the prefactors in 
Eqs. (13. 1611 and (13.171) . Hence we can normalize the two kernels such that the 
microscopic origin limit of these kernels is finite in preparation of section [5] We 
find the result for arbitrary 7 , 


det fe (X + tip) 
for k = 2m even and 

k 


V) / — C pk (t)pfo<a,6<fc-l 

P.7 


~(7.P+*0 

ab 


0 t ) 


det*(X + tip)) = C^(t) P f 0 < a , b < fc —1 

/ p ,7 


^^p+fc+l) 

‘ab 

_£(7,P+fc+l )( t ) 0 


(t) e 


( 7 ,p+fc+l) 


(t) 


for k = 2m + 1 odd. The functions serving as the kernels are 


Hi'V) = 


(-l)«+‘( 2 T+“+‘+ I l(! - 1)Z,*,,(<)„ 


Zi^(t) 


dZAfai - «2) 


(4.34) 

(4.35) 

(4.36) 


x (det (X - Kili_ 2 ) det (X - k 2 1 ;-2))J_2, 7 Ui=* 2 =-t 

_ ^_-^ja+h^ 27 +a+&+l 

(, v^ /2 d b K2 R^ ] («2,t) - («2, t) ( K l^) 


„Ct) / 


rr'(f) 


X 


and 


E 

i=o 

i e 2N, 

( (- 3)/ 2 ga^h) +i ^t) d b K2 R^ ] (k 2 ,*) - 3« 2 -R 2 }li («2,i)^“ 1 -R 2 ] ) («i,i) 


E 

J=0 


„(t) / 




l G 2N+ 1, 


(_1 'ioj- 27 +q 

#' 0 (f) = l ( /_ 2) , (det (X - «l i _ 2 ))J _ 2 ,7 U 


(4.37) 


i-1 i a +4 2 7+ a f ^-°-2^( _ U[7 + (l — l)/2,7 + I/ 2 , t/2] 

1 J \(l-a-2)\ U [7 + (Z — l)/2 ,7 + 3/2, t/2] 


^ 7 a-3 +1) (-t) 

(l — a — 3)! 


Due to this particular normalization we have the overall constant 


p! 


J p-\-k,~y 


>7 W 




(p+ky.n-Zof- tk(k+27)/2z p^ t y 


k G 2N 0 , 


P'- 


Z p -\-k— 1,7 (t) 


(4.38) 


j! t fc (fc+27)/2+(27-l)/2^ 7 (£)’ 


k G 2N 0 + 1. 


Note that we have not included the prefactors shown in Eqs. m and (ED for the 
full expressions of the gap probability and the distribution of the smallest eigenvalue. 
For the definition of the limiting quantities see Eqs. m and ED where all factors 
in p are accounted for. 

First we show the explicit expression of the gap probability at finite N. We 
multiply the results (14.341) and (14.351) for 7 = 0 with the factor e~ pt ^ 2 Z Pt o(t)/Z Pi 2 k- 
Then we find the first of our main results, 

TT u\ _ rWn„ 4 -\-k 2 / 2 + 1 / 2 r ,-pt/ 2 ^'[(p + k + 1)/2]U ((p+ k + l)/2,3/2; t/2) 

^pfikKP) — L 'p,fcl‘*W e 2y/2p> 


■' p,k 

X Pfo<a,b</c-l 


7j(0 ,p+k) 
ab 


(*) 


(4.39) 
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Figure 1. Visualization of our analytical results (14.39D , H4.40D , ( 14.42b , and (14.431) 
for finite p = 10 (solid curves) compared to Monte Carlo simulations (symbols, 
histogram) for the gap probability (left plot) and the distribution of the smallest 
eigenvalue (right plot). For the gap probability we generated 10000 real Wishart 
matrices of size 10 x (10 + z/) with an index v = 2k = 0,2,4,6, 8. For the 
distribution of the smallest eigenvalue we simulated 20000 random matrices with 
various dimensions p shown in the inset corresponding to v = 2k = 2,4, 6, 8. 


for even k and 


E pak (t) = Cp°l(4pt)~ k2 /2+1 e~ pt/2 F[(P + fc)/2]U fc)/2 ’ 3/2; t,2) 

£ v £p 


X pfd<a,b</c-l 


si° ,p+fc+1) (t) ^p+k+D^ 


_£-(0,p+/c+l) 

s b 


( t ) 0 


for odd fc, where we introduce the global normalization constant 


(4.40) 


(o) _ i - [' 4 i+1 (2/)!r[p + l + 2]p l 1 
dr[p + 21 +1] 

1 2" fc / 2 p!p 3fc/2 r[(p + l)/2] 
(p + fc)!r[(p + k + l)/2] 

2 -(fe+3)/2 p!p (ifc-i)/2 r ^ + 

. (p + fc)!r[(p + fc)/2] 


X J. 

1—0 


x 


, fc € 2N, 

, k € 2N + 1 


(4.41) 


which is also ideal to take the limit p —> oo. 

The explicit expression of the kernels E^ b ’ p+k ^ (t) and £(°>p+fc+i) ^ are no t m uch 
more enlightening at finite p than the expressions (14.361) and (14.371) for general 7 . 
Therefore we skip their expressions here and show them explicitly for the particular 
case 7 = 0 and 7 = 1 in the microscopic limit in section [5] For finite p we visualize 
the gap probability in Fig. |T] 

We skip the explicit kernels for the distribution of the smallest eigenvalue, too, 
and only explicitly show the Pfaffian structure with its normalization constant. For 
this reason we multiply the results (14.341) and (14.351) for 7 = 1 and p —> p — 1 with the 
factor pt^ 2k ~ 1 ' ) ^ 2 e~ pt ^ 2 Z p -i i i(t)/Z Pt 2 k from Eq. (13.91) 


P (f]dt _ C W ( ln , r fcV2 + i r -*/2 r I(P + k + 2)/2]U ((P + k + 2)/2,5/2; i/2) 

Jp,2fcfi;at — e 2(2 p) 3 / 2 


p,k 
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2 (y p+fe _ 1) (t )l d(4pt) 


(4.42) 
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for even k and 


PpM t)dt = F[(P + k + 1)/2]U 2( ( ^t/2 + 1)/2 ’ 5/2; ^ 


X Pfo<a,b<fc-l 


“ah \ l ) 

-d 1,p+fe) w 


d hp+k) (t) 

0 


cl(4pt) 


(4.43) 


for odd k. We multiplied the differential df to emphasize that this quantity is a density 
and transforms as a pseudo scalar under changes of coordinates. The constant 


k -1 


42=n 


i-i 


1=0 


A l+1 (2l)\T[p+l + 2]p l 

dr[p + 21 +1] 

1 2-( fe+ 5 )/ 2 p!p( 3 fc - 1 )/ 2 r[(p+ l)/2] 

(p + fe)!T[(p + k)/2] 

1 2"( fe+ 6 )/ 2 p!p 3 fc / 2 r[(p + l)/2] 


, k G 2N 
, k G 2N + 1 


(4.44) 


(p + A:)!r[(p + k + l)/ 2 ] 

it is t independent and converges to a finite number in the limit p —> oo. 

Equations (14.421) and (14.431) are our second main result. The distributions at finite 
p are visualized in Fig. [l] The two simplest cases for k = 0, 


P p , o(t) = 


p\ 


2P-i/ 2 T[p/2] ^ 


P pt/2 U ( 


P-1 

V 2 ’ 


1 t 
2 ’ 2 


(4.45) 


and fc = 1 


P p , 2(t) = s™U e -^ 


(4.46) 


U 


P-1 


1 t 

2 ’ 2 


(P- 1)! 


U 


+ 1 1 t\ (-l) p - 2 4 3 i 2 (-t) 


2 ’ 2 ’ 2 


(p — 2 )! 


obviously agree with the results by Edelman [25] . We underline that he has employed 
the standard normalization of the Laguerre polynomials while we have chosen the 
monic normalization. Moreover we have used Rummer’s identity of Tricomi’s confluent 
hypergeometric function and the duplication formula of the Gamma function to obtain 
Eqs. (14.45[) and (14.461) from Eqs. (14.421) and (14.431) . respectively. 


5. Microscopic Origin Limit 


In the microscopic origin limit of the gap probability (12.51) we have to perform the 
limit p —> oo while keeping 2k = v = n — p and t = u/4p fixed. Hence we zoom into 
a region of scale 1/p around the origin. In this region chiral random matrix theory is 
identical with physical theories like QCD below the critical temperature [Mi mi and 
condensed matter theory of disordered system for particular topological insulators and 
superconductors [5S]. In this regime the limiting gap probability and distribution of 
the first eigenvalue are defined as follows: 


and 


£ v (u) := lim E p v 

p—too 



(5.1) 


V u (u)du := lim P p „ 

p—too ’ 



(5.2) 
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We underline that the scaling factor resulting from the differential in the second 
definition is crucial to obtain a finite limit. 

The first question we have to address is the dependence of the limit on the parity 
of p (if p is even or odd) because the skew-orthogonal polynomials crucially depend 
on it, cf. Eqs. m and m- We can circumvent this problem by the fact that all 
important quantities can be written in terms of averages of characteristic polynomials 
to some powers, in particular the normalization (13.611 and the kernels (13.1611 and (13.1711 . 
Therefore it is sufficient to show that the limit of the average (ED is independent of the 
parity of p. For this reason we consider the intermediate result (14.711 . The microscopic 
limit can be readily performed by rescaling the integral U —> 2pU and r —> 2 pr yielding 

lim I p (—^ oc [ dfi(U)det J U [ drr( 27 ~ 2fe ~ 1) / 2 det 1 / 2 (rl 2 fe + 17) 
p->oo \4 pj J J Q 

CSE(2fc) 

x exp fi [tr Uk — tr f7 _1 + t r + r -1 ] j , (5.3) 


up to a p independent normalization constant. In particular this result is independent 
if the limit has been approached by an even p or an odd one. Therefore we do the 
asymptotic analysis for even p, only. 

The behaviour of the Laguerre polynomials and Tricomi’s confluent hypergeomet¬ 
ric function in the microscopic limit determines the whole asymptotics. The asymp¬ 
totics of both kinds of functions are given by 


and 


nr,- 1 M 

^o.r(«p + e+l) ip) 



b/2 


I b (y/au ), 


(5.4) 
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r(a f+ c) u 

p b ~l 
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, u 
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8 p 


= 2 



( b — 1)/2 

K b _i 



(5.5) 


respectively. Here we employ the modified Bessel functions of the first and second 
kind, I a and K a , respectively. The auxiliary parameters a,b,c take certain values in 
the particular limits below. 

The limits (15.41) and (15.51) imply the following asymptotics of the derivatives of 
the polynomials for even order 
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\C{xu) -2a- 2 ud u ] lim 

p-> oo \ 4 p/ 
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Here we have used the rescaled function (14.371) which is one of the polynomials for 
/c = 2m + 1 odd. The shift in the order of the polynomial from l — 2 to xp , cf. 
Eqs. (14.371) and (15.61) . has no effect on the limit since it can be absorbed in the limit 
p —> 00 . Furthermore we used the recurrence relation between the polynomials of 
even order and the polynomials of odd order, see Eq. (13.151) . to derive the limit for 
the latter. The ambiguous function c px / 2 {u/{4p)) is chosen such that its limit exists 
and converges to C(xu). Indeed this function only depends on the combination xu 
because a scaling in x yields an inverse scaling in u by absorbing the scaling factor 
in the p limit. We underline that also the constant C(xu) is ambiguous and we may 
choose it in such a way that the prefactor in front of the Bessel function I 2l + a +i(\fxu) 
vanishes, i.e. 
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(5.8) 


This expression is more compact than the one in Eq. (15.71) such that we stick with this 
intermediate result (which is again simpler than in [3J|). 

The normalization constants of the polynomials, see Eq. (14.151) . have the limit 


lim 


„(7) 


p ->00 T 2 [xp + l](a;p) 27 XP//2 \4p 


= 4. 


(5.9) 


Moreover we need to express the limit of the sum of p terms in order to deal 
with the two-point kernel (13.121) . Let us consider a function fj where the limit 
linip^oo f px = f(x) exists for all x £ [0,1]. Then the sum becomes 
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(5.10) 


Combining this limit with Eqs. (15.61) . (15.81) . 
kernel 
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Recall that the asymptotics for even and odd p yields the same answer, and so we 
have obtained this limit by choosing p = 2 L even. 

When considering the particular cases of the gap probability (7 = 0) and the 
distribution of the smallest eigenvalue (7 = 1 ) it is more enlightening to express the 
Bessel functions involved in terms of more explicit functions, 
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Then the gap probability reads 


/k -1 


£-2k(u) = ] 


4' + 1 (2/)! 


vZ =0 




-k 2 / 2-u/8-y774 


P f ( 


0<a,b<k— 1 


r-^(0,Oo) 

' 'ab 


(«) 


y/Ti 


pfo<o,b<fc-l 




(0,oo) 


(u) 


k G 2N 0 , 
k G 2N 0 + 1 


with 


= (j) a/ 2 [/o(^)+/o+i(^)] 
vV 2 


^ (0 j 00) 
~“ab 


(u) = 


da; 2 / 




2(b - a)I a (2x)I b (2x) 


(5.12) 

(5.13) 


(5.14) 

(5.15) 

(5.16) 


+ (26 + l)/ a+1 (2x)Jb(2x) - (2a + l)J a (2a:)/ 6+ i(2a:) 

The results in Ref. [52], have three typos. First, the factor x a+b+1 (^( a + b l / 2 in the 
notation therein) is missing in the integral for the two-point kernel see Eqs. (25) 
therein. Second, the indices of the Pfaffian in the case of odd k should go from 0 to 
k— 1, see Eqs. (27) and (29) therein. 

The gap probability seems to diverge at u —> 0 due to the terms u~ k / 2 and 
^-(fe 2 - 1 )/ 2 . However the kernels vanish as S^’°°^(u) oc u a+b+1 and oc u a for 

u < 1 such that both terms cancel and the gap probability behaves as a constant, 
especially it is normalized to f 2 fc( 0 ) = 1 . 

The distribution of the first eigenvalue is 
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Figure 2. The distribution of the smallest eigenvalue (|5.17[) (solid curves) 
smoothly fits with the microscopic level density (12.411 (dotted curves) around the 
origin. Only for larger values of u deviations due to the contribution of the second 
and larger eigenvalues become visible. A similar plot can be made for v = 0 which 
we have omitted for a better readability of the figure. 


The behaviour of the distribution at the origin can be read off from the kernels which 
are 


<«>=(?) 


/ 2+1 

y/u/2 
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Ia+2(Vu) H- -j= I a+3 (y/u) 
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I a+2 (2x)I b+3 (2x) 


Since (u) oc u a+b+3 and £i 1 ’°°' ) ( u ) oc u a+2 for |u| Clwe have T , 2k(.u) oc u k ~ 1 ^ 2 

agreeing with the behaviour of the microscopic level density (12.dl) . We remark that 
the term J u (y/u)/y/l&u ~ u^ u ~ 1 ^ 2 in the expression (12.41) is the dominant term in the 
limit u —> 0. The same behaviour of both distributions around the origin is inherent 
because the level density is governed by the smallest eigenvalue in this regime. Only 
for larger argument u the other eigenvalues start to contribute to the level density, cf. 
Fig. [2] 

The results (|5.141) and (15.171) are our third main result. We emphasize that these 
results do not only describe the smallest eigenvalue of an artificial system, namely 
real chiral Gaussian random matrices with an even index v = 2fc, but also of physical 
systems. Due to universality not only the level density (12.41) has to agree with those 
from physical systems like QCD or mesoscopic systems. Also the distributions of the 
smallest eigenvalues of those physical systems have to follow the same distributions of 
random matrix theory in the limit of its applicability, e.g. see [37(l38l (3] and references 
therein. This agreement should already happen at moderate system sizes. In Fig. [3] 
we compare the analytic results (14.421) and (14.431) of the distribution of the smallest 
eigenvalue at finite p with the microscopic limit (15.171) . This comparison underlines 
how fast the convergence to the universal result happens. Although this comparison 
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Figure 3. Illustration of the rate of convergence by comparing the finite p 
result 114.4211 for p = 11, 51,131 with the universal limit p —> oo 115.1711 for the 
distribution of the smallest eigenvalue. The index of the Wishart matrix is chosen 
v = 2k = 4 as an example. 


is done for random matrices which have particular advantages in comparison to real 
systems we expect that also physical systems should display a rapid convergence to 
the universal result. 

6. The Correlated Wishart-Laguerre Ensemble 

We are now interested in the effects of a fixed, non-trivial correlation matrix C not 
proportional to the identity matrix on the distribution of the smallest eigenvalue, see 
Eq. m- Such a correlation matrix can naturally encode system specific information. 
In time series analysis such a correlation may encode correlations between companies 
in finance m , seasonal effects in climate research [ 3 ] or organized crime in criminal 
defence [9]. However usually such correlation matrices have no microscopic limit in 
time series. But also in QCD and mesoscopic systems which exhibit such a microscopic 
limit correlations may appear. These correlations encode the structure of space-time 
and the choice of the gauge theory for example. They are system specific informations 
and may have an influence on the smallest eigenvalues. 

We expect that the correlation matrix has no influence on the smallest eigenvalues 
as long as its eigenvalues have a finite distance to the origin. Then the screening of the 
infinitely many eigenvalues between the smallest eigenvalue of W and the eigenvalues 
of C is strong enough. This was also shown in [l4j. For this purpose we choose 
a non-trivial empirical correlation matrix C / l p . In Fig. [3] we compare Monte- 
Carlo simulations with such an empirical correlation matrix C and the universal 
result (15.171) . The matrix size is chosen such that (n ~ p)/p = u/p -C 1 where 
p = 200 and v = 2k = 0,2,4. The perfect agreement underlines that correlations 
in the Wishart matrix have a very weak effect on the spectral statistics of the smallest 
eigenvalue. 

Again we emphasize that we have not looked at the situation where C develops 
a spectrum where some eigenvalues lie on the scale 1/p. Nor have we looked at the 
situation of doubly correlated Wishart-Laguerre ensembles. However for the latter 
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u 


Figure 4. Comparison of Monte Carlo simulations (symbols) with a non-trivial 
empirical correlation matrix C at finite p = 200, and the microscopic limit p — > oo 
of the distribution of the smallest eigenvalue 115.1711 (solid lines). The empirical 
correlation matrix is shown in the inset. We have generated 10000 correlated 
Wishart random matrices of sizes 200 X (200 + u) with v = 2k = 0, 2, 4. 

kind of ensembles we expect a similar if not exactly the same behaviour like one-sided 
correlated Wishart-Laguerre ensembles as considered here. 

7. Conclusion 

We have addressed and solved an open problem in the real Wishart-Laguerre ensemble 
also known as the chiral Gaussian orthogonal ensemble of rectangular p x n random 
matrices. We computed the distribution of the smallest eigenvalue and its integral, 
the gap probability that the vicinity of the origin is empty of eigenvalues. To this aim 
we have established that an integrable Pfaffian structure holds also when p — n = v 
measuring the rectangularity (or topology in the field theory application) is even. Such 
a Pfaffian structure was previously only known when v = 2k + 1 is odd. So far for 
an even rectangularity a recursive construction in p led to closed form expressions 
for v = 0,2 only. In view of the various applications it would be unnatural to 
restrict oneself to odd v, and not to expect for such an integrable structure to exist. 
However, the recursive construction (and closed results for v = 0,2) already revealed 
the appearance of special functions for finite p which are absent for v odd, namely 
Tricomi’s confluent hypergeometric functions. From our construction we now better 
understand why they appear through the expectation value of the square root of 
characteristic polynomials which are among the building blocks for the quantities in 
question. 

On a technical level our computation was possible due to the combination of the 
method of skew-orthogonal polynomials, though with a non standard weight, and the 
map of our building blocks onto invariant co-set integrals derived by bosonisation. 
In an initial step the computation of the the gap probability and the distribution of 
the smallest eigenvalue requires the evaluation of expectation values of characteristic 
polynomials raised to integer/half-integer powers for v odd/even, respectively. Hence 
the problem exhibits an increased level of difficulty in computing these objects for 
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even v. By including the square root for even v into the weight function we were back 
to expectation values of integer powers, which are known to be expressible through 
Pfaffian determinants of kernels and skew-orthogonal polynomials. The price we had 
to pay was to compute the latter for a non-standard weight including the square 
root. This was done by expressing the polynomials and kernel themselves through 
expectation values, mapping these back to matrix integrals and computing them via 
bosonisation. 

Indeed one can also consider the distribution and the cumulative distribution of 
the second to smallest eigenvalue, third to smallest eigenvalue etc. These quantities 
can be simply deduced from our results, too, because we calculated a quite explicit 
expression of the kernel at finite p. In particular we found a Christoffel-Darboux-like 
formula which expresses the sum over (p + k) terms in a sum over twelve terms, only. 

The Pfaffian structure enabled us to take the microscopic large-p limit at the 
origin, while keeping v = 2k fixed. In this limit we could show that the distinction 
between even and odd p for finite p becomes immaterial. We found results in terms 
of a Pfaffian comprising the limiting kernel for even k, plus an additional column and 
row for k odd, both for the gap probability and the smallest eigenvalue. 

Our results are universal for non-Gaussian potentials, as inherited from the 
universality of the known density correlation functions. We have checked that our 
findings follow the microscopic density for small argument, and that our fmite-p 
results, which we have confirmed through numerical simulations, converge towards 
the universal limit. Furthermore, we have also studied numerically the distribution 
of the smallest eigenvalue in an example of a correlated Wishart-Laguerre ensemble. 
We found that for moderate p it already follows the universal limiting distribution for 
several values of v. 

The computation of products of ratios of characteristic polynomials that also 
include square roots is in general an open question in random matrix theory. The 
structure of the results we obtained on a subset of such correlators should be relatively 
easy to translate to the Gaussian orthogonal ensemble, where such correlation 
functions enjoy further applications, e.g. in Quantum Chaos. It is very plausible 
that our universal result will also apply when introducing a fixed trace constraint, as 
it is known for odd v. What is less clear is whether a corresponding representation in 
terms of hypergeometric functions of matrix arguments exist, having the advantage 
that they can be continued to real /3 > 0. It would also be very interesting to see, 
if and when the universality at the origin breaks down for the correlated Wishart- 
Laguerre ensemble. A further open question is the computation of the gap probability 
and smallest eigenvalue distribution in the chiral Gaussian symplectic ensemble with 
/3 = 4. Apart from v = 0 only Taylor expansions exist so far. However, following 
similar ideas as in the present work, Pfaffian structures for these quantities exist and 
should be universal. Work in this direction is currently under way. 
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